Quantum Correlations in Two-Particle Anderson Localization 
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We predict the quantum correlations between non-interacting particles evolving simultaneously in 
a disordered medium. While the particle density follows the single-particle dynamics and exhibits 
Anderson localization, the two-particle correlation develops unique features that depend on the 
quantum statistics of the particles and their initial separation. On short time scales, the localization 
of one particle becomes dependent on whether the other particle is localized or not. On long time 
scales, the localized particles show oscillatory correlations within the localization length. These 
effects can be observed in Anderson localization of non-classical light and ultra-cold atoms. 



More than fifty years ago, P.W. Anderson predicted 
that a single quantum particle released in a disordered 
lattice can exhibit exponential localization in space [l| , a 
phenomenon termed Anderson Localization (AL). Since 
then there has been an ongoing effort to observe the sig- 
nature of AL experimentally, for example using light [2|. 
Recently a novel approach enabled the direct observation 
in space of AL for photons [3|-l6| and ultra-cold atoms 
0. These experiments, reporting the exponential local- 
ization of the particle density distribution, agree with 
the predictions of the single particle model as long as no 
interactions are involved. 

The question arises, whether there are measurable phe- 
nomena that are not described by the single particle 
model when considering the localization of a few pho- 
tons or atoms. It is known that when indistinguishable 
particles propagate together, exchange terms can result 
in the formation of correlations between their positions 
even in the absence of interactions. This result, known as 
the Hanbury Brown- Twiss (HBT) effect (I], was studied 
theoretically and experimentally for particles propagat- 
ing in free space. It was observed that quantum statis- 
tics play a crucial role in the formation of HBT correla- 
tions - bosons tend to bunch, while fermions exhibit anti- 
bunching ji|,[9[. Yet, it is not clear how HBT correlations 
evolve when the particles propagate through disordered 
media and become localized. While correlations were re- 
cently observed in the multiple-scattering of non-classical 
light experiments on AL focused on the parti- 

cle density distribution, and therefore could not observe 
spatial correlations between the localized particles. Now, 
experimental techniques that allow a direct observation 
of localization [1-0| can provide access to the study of 
HBT correlations between particles exhibiting AL. 

In this Letter we predict the correlations between two 
indistinguishable quantum particles evolving simultane- 
ously in disordered lattices. We consider non-interacting 
particles, and therefore the particle-density follows the 
single particle dynamics: both particles undergo Ander- 
son localization. Nevertheless, we find that the two par- 
ticles develop non-trivial spatial correlations due to in- 
terferences of all the scattering paths the two particles 
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Figure 1: (color online). Dynamics of a single quantum parti- 
cle placed at t = on a single site, (a) In a periodic lattice, the 
particle density distribution expands ballistically, with high 
density at the edges of the distribution, (b) Cross-section of 
the density distribution after some evolution. (c),(d) In a dis- 
ordered lattice, the expansion is limited to a finite region, (e) 
The density distribution, averaged over 1000 realizations of 
disorder. The expansion starts ballistically. After some prop- 
agation a localized component emerges around the initial po- 
sition, (red arrow, red cross-section in (f)). The ballistic com- 
ponent decays in time, leaving the distribution exponentially 
localized (blue arrow, blue cross-section in (f)). 



can take as a pair. On short time scales, the localiza- 
tion of one of the particles uniquely determines whether 
the other particle will be localized or not. On longer 
time scales, when both particles localize, the particles 
exhibit oscillatory correlations within the localization 
length. Remarkably, these oscillatory correlations sur- 
vive multiple scattering even after very long evolution 
times. In addition, we show that fermionic correlations 
can be reproduced by entangled bosonic states. 

To study two-particle localization, we follow the for- 
malism developed in [13| for periodic systems. We con- 
sider a one dimensional quantum tight-binding model, 
given by Hamiltonian: 
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where a^ n is the creation operator for a particle in site n, 
W n is the on-site energy and T n?m is the tunneling ampli- 
tude between nearest neighbors. At time t the creation 
operator at site r is given by 



(a) Bosons 



ai(t) = J2Urr>(t)al,(t = 0) 
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where the time-evolution operator U rr >(t) = (e' lHt ) r ^ is a 
unitary transformation given by calculating the exponent 
of the Hamiltonian matrix H, and describes the ampli- 
tude for transition of a single particle located at site r' 
at t = 0, to site r at time t. In the following, we study 
the evolution of two indistinguishable particles in a dis- 
ordered lattice, where disorder is introduced by random- 
izing the tunneling amplitudes. We focus on the particle- 
density n r (i)=(ala r \ and on the two-particle correlation 
Fq jr (t) = (a^ q ala r a q ), highlighting the single-particle ver- 
sus two-particle features of the dynamics. 

We start by considering the evolution of a single parti- 
cle, described by the particle density n r (t) = \U rr >(t)\ 2 , 
where r' is the location of the particle at t = 0. In Fig. 1 
we plot the evolution of the particle density for an initial 
condition aj |0) , i.e. a single particle at site number 0. 
Fig. la depicts the propagation of the particle in a pe- 
riodic lattice, which is given by n r = J^_ r ,(2Tt), where 
J n is the Bessel function of order n. This dynamics is 
known as a quantum walk (l3l-[l6|. in which the particle 
density expands ballistically (i.e. linearly in time), and 
the probability to find the particle is highest at the edges 
of the distribution (see Fig lb). When disorder is intro- 
duced (here by randomizing the tunneling amplitudes), 
the propagation changes significantly, as shown in Fig. 
lc,d. A statistical average of the probability distribu- 
tion over many realizations of disorder (Fig le,f) reveals 
a clear cross-over from ballistic expansion at short times 
to exponential localization at long times (blue arrow in 
Fig le). At intermediate times (red arrow in Fig. le), 
the particle can be either localized or ballistic 

We now consider two non-interacting particles which 
are initially localized on sites r' and q' focusing on the 
correlations between their locations at a later time t. 
The correlation matrix is given by Y q ^ r = {{\U qr >U rq > ± 
U qq 'U rr '\ 2 )) (H)]. Here ((•)) signifies averaging over re- 
alizations of disorder, and the upper and lower signs de- 
scribe, respectively, bosons and fermions of same spin. 
The depicted matrices (for example in Fig. [2j) represent 
the probability to detect at time t exactly one particle 
at site q and one particle at site r. As we show below, 
the density distribution is identical for bosons and for 
fermions (see Fig. [2^), yet the emerging quantum corre- 
lations between the particle pair depends on the initial 
position of each particle and on their quantum statistics. 

We first consider the case in which initially the two par- 
ticles are positioned at two adjacent sites in a separable 
state, i.e. the initial condition is given by \(j>) = aja^O). 
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Figure 2: (color online). Two-particle correlation function in 
disordered lattices. At t = 0, two particles are placed at two 
neighboring sites. The correlation function is calculated after 
a short evolution, where remnants of the ballistic component 
are still present (see Fig. le,f), averaged over a 1000 disorder 
realizations. The matrices represent the probability to detect 
one particle at site q and one particle at site r. (a) Bosons 
either scatter and become localized, or remain ballistic. How- 
ever, if they both remain ballistic, they always reach the same 
side of the distribution, (b) Magnification of the localized 
component at the center of the correlation map. (c) Fermions 
show similar behavior, only that if they both remain ballistic 
they always separate to different sides of the distribution, (d) 
Magnification of the correlation map, showing a checker- like 
pattern, (e) The particle density distribution, which is identi- 
cal for fermions and for bosons, computed by tracing out the 
position of one of the particles (summing over the columns of 
the correlation matrix), (f) The inter-particle distance prob- 
ability for bosons (blue) and fermions (red). Bosons tend to 
appear at the same site, while fermions are more likely to be 
separated by an odd number of sites. 



We calculate the two-particle correlation function after 
a relatively short evolution time, in which each particle 
has a non-zero probability to be localized or to remain 
ballistic (see Fig. 1). As shown by the correlation ma- 
trix in Fig. 2a, if the two particles are bosons they can 
both remain ballistic (corners of the correlation matrix), 
or both become localized (center of the correlation ma- 
trix). The possibility for one boson to localize and the 
other to remain ballistic is also non-zero. However, when 
both bosons exhibit ballistic behavior, the probability 
for the two bosons to separate to two different sides of 
the distribution vanishes (top left and bottom right cor- 
ners of the correlation matrix) . This is the signature of 
bosonic bunching in lattices [13] . In this case then, the 
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Figure 3: (color online). The two-particle correlation func- 
tion when At £ = 0, the two particles are placed at two non- 
adjacent sites, (a) Bosons show two possible effects - they 
either both localized, or both remain ballistic, (b) Closeup on 
the correlation inside the localized component, (c) Fermions 
almost always split - one fermion becomes localized and the 
other remains ballistic, (d) Closeup on the correlation inside 
the localized component, (e) The particle density distribu- 
tion, identical for fermions and bosons, (f) The inter-particle 
distance probability within the localized regime for bosons 
(blue line) and fermions (red line). The fermions exhibit a 
flat inter-particle distance distribution (except at zero sepa- 
ration), while the bosons show an oscillating probability. 



two-particle correlation has two components - a ballis- 
tic component showing spatial bunching, and a localized 
component without spatial correlations. A closeup on 
the correlation in the localized component is given in Fig. 
[2]3, also describing the correlation matrix at later times, 
when the distribution is completely localized jl7j [. 

In the case of fermions, the correlation matrix of Fig. 
[2k, d shows that if the two particles remain ballistic (i.e. 
both particles did not localize yet), they will exhibit anti- 
bunching and separate to opposite sides of the distri- 
bution. A closeup of the correlations inside the local- 
ized component reveals a checkers-like pattern, mean- 
ing that even when both fermions are localized they are 
non-trivially arranged in space. Most significantly, this 
checker-like pattern survives the random scattering pro- 
cess even after very long evolution times jl7j |. 

To understand the significance of the checker-like pat- 
tern we extracted from the correlation matrix the inter- 
particle distance probability given by g(A) = r^+A- 
Comparison between the results for bosons and fermions 
are depicted in Fig. [2)P. Bosons (blue line) have the high- 



est probability to be at the same site, and the probability 
drops monotonously with distance. In contrast, fermions 
have zero probability to be on the same site (as expected 
due to Pauli exclusion). Yet, interestingly, the correla- 
tion oscillates: the two fermions tend to be separated by 
an odd number of sites. 

An even more dramatic effect takes place when the ini- 
tial condition is such that the two particles are placed 
at no n- adjacent sites. Consider the initial condition 
\4>) = a}_ l a\ |0) - two particles at two sites separated by an 
empty site. Again we calculate the correlation function 
after a relatively short propagation, in which the dynam- 
ics includes both a ballistic and a localized component. 
For bosons (Fig. [3^) there are only two possibilities: ei- 
ther the particles are both localized (center of the corre- 
lation matrix) OR they both remain freely propagating 
(corners of the correlation matrix). The possibility for 
one particle to be localized and the other one to be free 
is diminished to nearly zero. A closeup on the correlation 
when both particles are localized (Fig. [SJd) shows that in 
this case it is the bosons that exhibit the checkers-like 
pattern; they are more likely to be separated by an even 
number of sites. 

For this input state, fermions exhibit a different sym- 
metry in the correlations (Fig. [St): the two particles tend 
to always occupy different components of the expanding 
wavepacket: one fermion becomes localized, while the 
other fermion remains freely propagating in the ballistic 
component. The closeup on the localized regime (Fig. 
[3]i) shows a flat distribution, except for the zero proba- 
bility to be at the same site. The inter-particle distance 
for this case is depicted in (Fig. [3F) [17] . 

Finally, it is of interest to consider the propagation 
in disordered lattices of path-entangled input bosonic 
states, which can be regularly generated in quantum- 
optics experiments. In a lattice system, such a state is 
given by a superposition of two bosons - placed together 
at one site OR another site \<j>) = ^[(aj,) 2 + e^(a£,) 2 ]|0). 
Interestingly, properly constructed entangled states re- 
produce some of the features presented by fermions. In 
Fig 0] we present the calculated correlation matrices for 
two kinds of path-entangled input states: In Fig. Hk, 
we plot the correlation matrix when the initial state is 
\<p) = ^[(«o) 2 + ( a l) 2 ]|0)- The correlation matrix shows 
several features that are similar to the ones exhibited 
by the fermions as seen in Fig. : the two parti- 
cles can both become localized, both remain ballistic, 
or split - one localized and one ballistic. However, they 
never end up on the same ballistic lobe. The main dif- 
ference is in the localized region: the checker pattern 
disappears, and the diagonal q = r of the correlation 
matrix now shows enhanced probability. In Fig. [4fc 
we plot the correlation matrix when the initial state is 

10) = ll( a -i) 2 + ei7r ( a i) 2 ]l°)- i- e - the two states are 
separated by an empty site. This time the pair shows 

a tendency to split - one boson is localized, while the 
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Figure 4: (color online) Correlation maps for path-entangled 
bosons evolving in disordered lattices, (a) Correlation map for 
the initial condition \4>) = |[(aj) 2 + (a\) 2 ]\0). (b) Magnified 
view of the localized component, (c) For = ^(al-J 2 — 
(aj) 2 ]|0). (d) Magnified view. These correlations are similar 
to the fermionic results in Fig. 2c, d and 3c, d (see text). 



other remains ballistic. The correlation matrix is similar 
to that exhibited by the fermion pair in Fig. [3b, except 
the enhanced probability on the matrix diagonal q = r. 

The emergence of the checker-like patterns in the cor- 
relations can be explained in terms of the lattice eigen- 
modes excited in each case. Disordered lattices support 
two kinds of eigenmodes: flat-phased and staggered, in 
which adjacent sites are in phase /tt out of phase, corre- 
spondingly jij. Certain initial conditions, (depending on 
the initial distance between the two particles and their 
quantum statistics) involve the simultaneous excitations 
of both kinds of eigenmodes, resulting in a density pat- 
tern that contains a component with two-site periodic- 
ity. This effect is washed out in the density distribution 
averaged over all realizations of disorder, as in each re- 
alization the oscillations appear in a different location. 
However, the fact that such oscillations appears in each 
realization will be recorded in the averaged correlation. 

At higher dimensional disordered systems, the single 
particle dynamics can be qualitatively different. For ex- 
ample expanding wavepacket can exhibit diffusive expan- 
sion before localization [5]. Nevertheless, we have verified 
that the results described above appear also in two di- 
mensional lattices, including the checker-like patterns. 

In conclusion, we studied the evolution of two non- 
interacting quantum particles in disordered lattices ex- 
hibiting Anderson localization. The two-particle corre- 
lation exhibits unique features that survive the random 
scattering, such as oscillating correlations and correlated 
cross-over from ballistic propagation to localization. The 
correlations depend on the initial separation between the 



particles and on their quantum statistics. Experimen- 
tally, the bosonic correlations can be observed by inject- 
ing photon pairs generated using spontaneous paramet- 
ric down conversion, into disordered waveguide lattices. 
Such lattices were recently used to observe Anderson lo- 
calization of classical light 0, 0| , and were shown to be 
feasible for observin g no n-classical correlations between 
photon pairs [HI, [H, H)]. With atom-matter waves, cor- 
relations can be observed in the quantum walks [16| of 
a pair of atoms in disordered lattices, offering a possibil- 
ity to study the effect of interaction [20| on correlations. 
Alternatively, Density-density correlations can be mea- 
sured for expanding Bose-Einstein condensates in disor- 
dered optical potentials 0, EH • 
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